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We report first principles characterization of the structural and electronic properties of 
(SrTi03)s/(SrRu03)i superlattices. We show that the system exhibits a spin-polarized two- 
dimensional electron gas, extremely confined to the Ad orbitals of Ru in the SrRuOs layer. Every 
interface in the superlattice behaves as a minority-spin half-metal ferromagnet, with a magnetic 
moment of n — 2.0 £tB/SrRu03 unit. The shape of the electronic density of states, half metallicity 
and magnetism are explained in terms of a simplified tight-binding model, considering only the t^ g 
orbitals plus (i) the bi-dimensionality of the system, and (ii) strong electron correlations. 



PACS numbers: 74.70.Pq, 73.20.At, 75.70.-i, 75.10.-b 

The family of ABO3 perovskite compounds is really 
prominent among all the complex oxides. The delicate 
interaction between their electronic, spin, lattice and or- 
bital ordering degrees of freedom, whose respective en- 
ergy couplings arc all of the same order of magnitude pQ , 
gives rise to a wide variety of ground states and phenom- 
ena [2]. These compounds typically exhibit exceptional 
functional properties that, moreover, can be finely tuned 
by playing with temperature and/or electrical (screen- 
ing) [3] , mechanical (strain) [3] , and chemical (doping) [5] 
boundary conditions. 

At the bulk level, ABO3 compounds cover the whole 
spectrum of conducting properties, ranging from good 
insulators and semiconductors to metals and even super- 
conductors ,2j. Moreover, it has recently been shown that 
their conducting behaviour can also drastically change at 
surfaces [B] and interfaces [7]. The possibility to create 
highly-confined two-dimensional electron gases (2DEGs) 
at oxide interfaces has also been demonstrated and gen- 
erated a huge excitement. 

The most widely studied system is certainly the 2DEG 
formed at the LaO /Ti02 polar interface between LaA103 
(LAO) and SrTi0 3 (STO) [8], two good band insulators 
at the bulk level. Although the mechanism at the origin 
of the transfer of charges remains partly controversial, 
the 2DEG is formed by the accumulation of electrons 
in the Ti 3d states of STO over a thickness of a few 
nanometers near the interface, and exhibits interesting 
exotic properties. In a different context, when looking 
for oxides with enhanced thermoelectric properties, Ohta 
et al. [9 highlighted the possibility to create a highly- 
confined 2DEG in [(STO) m /(SrTi . 8 Nbo.20 3 )„]p super- 
lattices. Since Nb-doped STO is a well-known conducting 
electrode material |10j , such a superlattice can be viewed 
as the limiting case of ultrathin metal oxide layers em- 
bedded in an insulating STO matrix. The 2DEG arises 



from the doping produced by the Nb atoms so, alter- 
natively, the system can also be interpreted as a partial 
doping of a STO matrix at the B site. Very recently, 
the formation of a 2DEG has been achieved through a 
related route: the doping of the A site of the STO ma- 
trix. The Sr atoms of a single SrO layer of STO are 
partially [TT] or completely [T^] replaced by a rare-earth 
element, R, which provides additional electrons. In this 
case, the conduction electrons provided by the substitu- 
tional layer are transferred to the STO matrix but stay 
near the RO layer due to Coulomb attractions, and pos- 
sibly also to the large dielectric constant of STO [T3] . 

In this Letter, we predict from first principles cal- 
culations the appearance of a highly-confined 2DEG 
in (STO) 5 /(SrRu0 3 )i superlattices (STO/SRO). This 
system is to some extent related to those previ- 
ously considered. SRO is a well-known conductive 
material so that STO/SRO present similarities with 
(STO) m /(SrTi . 8 Nb .2O 3 )„ superlattices 0. Also, 
STO/SRO can be seen as a STO matrix in which a Ti02 
layer is periodically replaced by a single Ru02 layer, ap- 
pearing as a variant of the work of Jang et al. [12] . How- 
ever, contrary to all recent examples of 2DEGs in per- 
ovskites previously discussed, that confine the conduct- 
ing electrons in Ti 3d orbitals of STO, the conducting 
electrons in STO/SRO are localized exclusively in Ru 
Ad orbitals. We will show that this gives rise to new 
and unexpected features. Contrary to current though 
regarding SRO thin films [14H16J . and finite period su- 
perlattice [17] . our infinite period superlattices exhibits 
an extremely confined half-metallic ferromagnetic (FM) 
2DEG, very promising for spintronic devices [18]. 

For this study, we perform first principles simulations 
of STO/SRO superlattices using two complementary ap- 
proaches. Unless otherwise stated, all the calculations 
presented here have been performed within the local spin 
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density approximation (LSDA) to the density functional 
theory using the Siesta code [19] . An extra Hubbard-U 
term is included to account for the strong electron cor- 
relations, with a {/ ff of 4.0 eV applied only to the d 
orbitals of Ru. The robutness of the results presented 
below have been doubled checked using the Crystal09 
code [SUJ within the Bl-WC hybrid scheme [ST]. Hy- 
brid functionals are well known to provide improved de- 
scription of magnetic and highly-correlated oxides, and 
Bl-WC appeared recently as a successful approach to 
describe 2DEG at the LAO/STO interface [22 . The 
excellent agreement between the results obtained inde- 
pendently within both approaches strongly supports our 
predictions. The procedure to fit the U e g and other tech- 
nicalities can be found in the Discussion 1 and 2 of the 
Supplemental materials |23| . Both spin-orderings in the 
SRO layers, FM and antiferromagnetic (AFM), where 
the up and down spins of the RuC>2 layer are ordered in 
a checkerboard arrangement, are simulated. 
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FIG. 1. (color online) Schematic representation of the re- 
laxed structure of STO/SRO superlattices with LDA+U. (a) 
Unit cell periodically repeated in space. O atoms are shown in 
red, Sr in green, Ti in blue and Ru in grey, with O octahedra 
colored accordingly, (b) Rotation angle of the oxygen octahe- 
dra, 9, around the [001] direction (z axis). Squares represent 9 
for the corresponding octahedra at the same height as in panel 
(a). Dotted lines are the theoretical value of 9 for bulk STO. 
(c) Layer-by-layer rumpling, r\i = [z(M;) — z(0;)]/2, where 
z(Mi) is the position of the cation, and z(Oi) is the position 
of O, in a given layer i. Atoms in layer i move in the direc- 
tion indicated by the arrows. Black lines represent the mean 
position of each atomic layer. Numbers inside the structure 
are the change in the interplanar distance between consecu- 
tive planes with respect to the ideal unrelaxed structure (half 
of the bulk value of STO, a — 3.874 A). Similar results are 
obtained within the Bl-WC fuctional (see Discussion 3 of the 
Supplemental material [23] )• Units in A. 



In the most stable phase, that turns out to have FM 
order, our simulations predict a PA/mbm space group. 



Oxygen octahedra rotate changing the phase from layer 
to layer, as shown in Fig. |T](b) . This pattern of rotation 
shares some similitudes with the a°a°c~ Glazer rotation 
of the STO matrix, but now the module of the rotation 
angles around the z-axis are not the same from one oc- 
tahedra to the adjacent one. The angle of rotation of the 
Ru06 octahedra, #sro = 11-6°, is roughly the same as 
in bulk tetragonal SRO under the same strain condition 
in the basal plane (0|gg = 11.4°). The rotations of the 
TiOg octahedra rapidly converge to the bulk LDA value 
of STO (6^ = 4.3°), indicating that the rotations are 
functions of the local chemical environment [24] . Tak- 
ing into account that bulk SRO crystallizes within the 
GdFe03 structure (Puma space group), with a~b + a~ 
Glazer rotation, we have explored the possibility of the 
existence of out-of-plane tiltings in SRO. However, un- 
der the epitaxial configuration considered here, we have 
checked that they dissappear. This supports that the 
single unit cell of SRO in our superlattices are stabilized 
in a pseudocubic perovskite structure |25j . As shown in 
Fig. [ljc), both the rumplings, and the changes in inter- 
planar spacing with respect to the ideal (non-relaxed) 
structure, are very small. 

Figure [2] shows the projected density of states (PDOS) 
on the atoms at the different SrB03 unit cells (where B 
is Ti or Ru), within the LDA+U [panel (a)], and Bl- 
WC hybrid functional [panel (b)]. The physical con- 
clusions that can be drawn are robust, independently 
of the theoretical framework used. In particular: (i) 
remarkably the system presents a 2DEG, strictly con- 
fined in the SRO layer. The electronegativity of Ru 
(2.20), larger than the one for Ti (1.54), is the main 
candidate to be the driving force for the electron local- 
ization. This localization resembles that recently found 
on the d orbitals of V in SrV03/Nb:SrTi03 quantum 
wells [26]. (ii) The electronic states around the Fermi 
level are fully spin-polarized, with only the minority- 
spin electrons involved in the charge transport |18j . (hi) 
The 2DEG exhibits ferromagnetism, with a magnetic mo- 
ment of \i = 2.0 /ie/SRO formula unit. Integrating the 
spin-polarized charge density on spheres surrounding the 
atoms (in Siesta) or analyzing the Mulliken spin popu- 
lation within Bl-WC (in Crystal), we see that the mag- 
netic moment is mostly due to the Ru atoms, /l/,r u = 1.4 
/ib, within the range of measured magnetic moments in 
single crystals [ST] , and polycrystalline samples 28 . The 
same value of /zr u is obtained in the AFM configuration. 
The remaining magnetic moment in the FM configura- 
tion comes from the oxygen atoms. These results are 
robust with respect to pseudopotentials, basis set, func- 
tional, and inclusion of spin-orbit coupling (Discussion 
4 of the Supplemental material 23J ) . Interestingly, the 
electronic structure at the interface inherits some of the 
particularities theoretically predicted in bulk SRO when 
the electronic correlations are properly included [as in 
LDA+U JSpJ, pseudo self- interaction correction [3D], or 
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FIG. 2. (color online) Layer-by-layer PDOS on the atoms at the SrBC>3 layers (B = Ti or Ru), within (a) the LSDA+U, 
and (b) the Bl-WC hybrid functional, for the corresponding relaxed STO/SRO superlattice. Each panel represents the PDOS 
of a unit cell, numbered as in Fig. [ija). Corresponding layers in the upper-half of the structure are equivalent by symmetry. 
Majority/minority spin are represented in the upper/lower halves of each panel, respectively. 



Bl-WC hybrid functional (Discussion 5 in the Supple- 
mental materials [13]). Between these properties, we 
might cite the half-metallicity [25] or the saturation of 
the magnetic moment [30]. Experimental observation of 
the half-metallicity remains a challenge, mostly due to 
the large magnetic fields required to overcome the high 
magnetic anistropy barriers and fully magnetize the sam- 
ples. 

Regarding the band shape, we observe in Fig. [3] that 
the bands with Ku-d xy character are significantly dif- 
ferent from the degenerate Ru-d XZyVZ bands. While the 
PDOS of Ru-d xy resembles that usually found in bulk 
t2 g bands [31] . the double peaked shape [peaks marked 
with arrows in Fig. gb)] for Ru-d xz,yz is similar to that 
of a ID metal. It is also significant that the relative po- 
sition of these bands for the majority and minority spins 
is different, an effect attributable to the strong electronic 
correlations in the system, enhanced with respect to bulk 
due to the two-dimensionality of the system (see below) . 
Indeed, to study the effect of the Hubbard-U in our calcu- 
lations we carried out a plain LSDA simulation (without 
U corrections) and plotted the PDOS (Discussion 2 in 
the Supplemental materials 23J). The shape of the Ru- 
d xy and Jhi-d X z,yz is the same as when the Hubbard-U 
term is considered. However, the position of the bands 
is shifted to almost the same energy range for both spin 
channels, destroying the half-metallicity. 

If the half-metallicity can be attributed to the elec- 
tronic correlation, the shape of the bands can only be 
explained through the 2D character of the electron gas. 
To check this point, we have made use of a tight-binding 
model, where only the ti g states are retained in the basis 
set. Qualitatively, the shape of the t2 g bands is found 
to be the same for all perovskites including transition 
metal ions in the B-position. Assuming cubic symme- 
try, the PDOS associated to these bands for a full three- 
dimensional solid is shown in Fig. pta) . Using this model 



we then simulate the changes in the PDOS when the sys- 
tem is bidimensional, extended periodically through the 
xy plane. In Fig.|4]jb) we can see that, while d xy -band re- 
mains unaltered, the d xz>yz ones now present two-peaks 
in good agreement with the ones in the full first prin- 
ciples calculations [arrows in Fig. [3jb)] . In particular, 
the PDOS of bidimensional d XZjyz levels closely resem- 
bles that of textbook one-dimensional electron gas [3T] 
due to the negligible interactions of d xz ^ yz wavefunctions 
in the xy plane with other than first neighbors. However, 
we can observe that the LDA+U bands are shifted with 
respect to the ones obtained in our tight-binding model. 
These shifts are due to (i) the difference in bonding of the 
orbitals with in-plane and out-of-plane neighbors and (ii) 
the Hubbard-U correction. To take into account the first 
effect we note that the d xy wavefunctions only interact 
strongly with in-plane neighbors which are all functions 
centered around Ru 4+ ions. In contrast, d xz ^ yz wavefunc- 
tions interact in-plane with centers around Ru 4+ ions 
while out-of-plane they interact with functions centered 
on Ti 4+ ions mediated by O atoms [see Fig. Be)]. Since 
Ru 4+ ions are more electronegative than Ti^ ones, the 
apical oxygens are polarized towards the Ru 4+ . This 
translates into an increase of covalency and, therefore, a 
decrease in the energy of the bonding states and an in- 
crease of the antibonding levels, denoted as A_E cov in Fig. 
Qc). Finally, we include the effect of the U e g value to 
take into account magnetism. LDA+U theory predicts 
that orbital energies are shifted by U c s (1/2 — A) where A 
is the occupation of the orbital. Using this formula, and 
assuming that the majority spin ti g levels are full, while 
the minority spin d XZl d yz and d xy contain respectively 
0.5, 0.5 and 0.0 electrons, we find the shifts depicted in 
Fig. Bd). Qualitatively, we find the bands to be in good 
agreement with those obtained in the full first-principles 
calculations [see Fig. |3jb)]. 

In summary, our first principles simulations pre- 
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FIG. 3. (Color online) PDOS on the atoms at the central 
SRO layer, projected on the orbitals of (a) an apical O atom, 
(b) the Ru atom, and (c) an equatorial O atom. In the DOS 
projections corresponding to Oxygen, the 2p orbital directed 
towards the transition metal ions is denoted by O(a), while 
the two quasi-degenerate orbitals perpendicular to this bond 
direction are O(tt). The arrows in panel (b) mark the position 
of the peaks for the Ru-d XZtyz orbitals (see text). 




FIG. 4. (Color online) Result of the tight-binding model used 
to interpret the DOS of i2 9 d orbitals. Our model shows 
the transformation from (a) the three-dimensional bulk, to 
(b) a two-dimensional structure. Refinements of the model 
to include, first, covalent bonding effects in-plane and out- 
of- plane (c), and second, shifts of the bands due to strong 
electronic correlation via a Hubbard-U term (d), are required 
to reproduce the first principles band structure. Black lines 
in (a) represent the d xy<XZtyz manifold, degenerated in energy. 
Once the symmetry is broken, the d xy state is represented by 
a green line, and the degenerated d X z,yz orbitals by red lines, 
(e) illustrates the orbitals used in the tight-binding model 
and differences of in-plane and out-of-plane bonding for d xz 
orbitals (identical sketch can be done for the d yz orbitals). 



diet the appearance of a two-dimensional, half- 
metallic, ferromagnetically ordered electron gas 
at the insulator /ultrathin- metal-film interfaces in 
(STO) 5 /(SRO)i superlattices. The 2DEG is extremely 
confined at Ru Ad orbitals. At odds with previous 
realizations of 2DEG, where the main mechanism 
behind localization is of electrostatic nature [T31 [55] , we 
propose the larger electronegativity of Ru to be the main 
cause for the extreme confinement. The half-metallic 
properties and the concomitant FM is an intrinsic 
property due to enhanced electron correlation, and not 
to extrinsic properties (O vacancies) as postulated in 
LAO/STO 2DEG 32J. Our results encourage exper- 
imental verification, to check the usefulness of these 
superlattices in magnetoresistance or spintronic devices. 
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